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ABSTRACT 


_ The total nuclear reaction cross sections for 180 MeV protons incident on Li, Be, C, Al, Ca, Fe, 
In and Au have been measured. Range separation was used to distinguish between elastically 
and inelastically scattered paiticles and the effect of the shape of the range curve was studied. 
The total rms errors are 2-3 % and mainly due to counting statistics. 

The experimental cross sections have been compared with the results of a calculation using an 
optical model with parameters derived from the Gammel-Thaler nucleon-nucleon scattering ma- 
trix. 


- 1. Introduction 


The optical model which has been commonly used to analyse experimental data 
on the elastic scattering of protons from nuclei is also capable of predicting the total 
reaction cross section. In this connection a “‘reaction’’ means an interaction in which 
the proton leaves the nucleus in an excited state or transforms it into another nucleus; 
i.e. all nuclear interactions are included except elastic scattering. The measurements 

‘reported in this paper were made at a bombarding energy of 180 MeV with the main 
purpose of reducing the ambiguity in the optical model analysis of elastic scattering 
at that energy. In particular, one would expect the imaginary part only of the 
potential to determine the reaction cross section quite well, whereas all terms con- 
tribute to the differential cross sections and polarizations. The details of the analysis 
are reported in the following paper [1]. 

Total reaction cross sections for nucleons incident on nuclei have been studied 
previously and the experiments concerned with protons as bombarding particles at 
energies between 30 and 900 MeV are listed in Table 1 together with some relevant 
experimental parameters. In all cases the attenuation in intensity of a proton beam 
by an attenuator (target) was measured, but the experiments may be divided into 
two groups according to the different methods of detecting protons emerging from 
the target. In the 134, 860 and 895 MeV experiments a detector that was insensitive 
to the proton energy was used and the reaction cross sections were obtained by extra- 
polation to zero solid angle of the poor geometry loss cross sections. Certain assump- 
tions regarding the angular distributions of elastically and inelastically scattered 
particles had to be satisfied, which made the extrapolation procedure somewhat 
uncertain and the method not very easy to use when aiming at accuracies of the 
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Table 1. Measurements of the total reaction cross section for protons at bombardin, 
energies between 30 and 900 MeV. 
a oe a ie eR 


Level below 


Mean proton | Energy band| which inelastic Elements esadied Ref. 
energy (MeV)| width (MeV) events are 
excluded (MeV) 


(a a 


~~, 


sme DP Yep 


34 5 0 C, Al, Fe, Sn, Pb Gooding [2] 

61 5 0 C, Al, Fe, Sn, Pb Meyer et al. [3] 

Ga 10 0 | C, Al, Cu, Ag, Pb Goloskie, Strauch [4] 

95 10 0 C, Al, Cu, Ag, Pb Goloskie, Strauch [4] 
113 10 0 C, Al, Cu, Ag, Pb Goloskie, Strauch[4]| 
133 10 0 C, Al, Cu, Ag, Pb Goloskie, Strauch[4] 
134 26 0 C, Al, Ga, Cd, Pb Cassels, Lawson [5] 
180 10 0 Li, Be, C, Al, Ca, Fe, In, Au] This experiment 
185 45 20 Be, C, Al, Cu, Pb, U Kirschbaum [6] 
240 65 20 Be, C, Al, Cu, Pb, U Kirschbaum [6] 
305 65 20 Be, 'G, Al Ca; Pb, Ui Kirschbaum [6] 
230 220 15 U Millburn et al. [7] 
290 100 15 Cy Al Cus Millburn et al. [7] 
860 60 0 Be, C, Al, Cu, Sn, Pb Chen et al. [8] 
895 70 0 C, Al, Cu, Sb, Pb Booth e¢ al. [9] 


order of a few percent. In principle, however, the separation of elastic and inelastic 
events was complete. The use of an energy-sensitive detector for the emerging parti- 
cles—as in all other experiments listed—introduced a means to set a separation level 
for inelastic events, i.e. to discriminate against protons having lost a certain amount 
of energy in the target in excess of the ionization loss. Thus the reaction cross sections 
could be obtained more directly. The separation level is indicated in the table, 0 MeV 
meaning that the inelastic contribution has been corrected for. As the protons lose 
energy in the target by ionization the value obtained for the cross section is an 
average value over an often considerable energy range. This may not be too serious, 
however, as the energy dependence of the cross sections seems to be slight. Table 1 
shows that when the present experiment was started there existed an early deter- 
mination of the reaction cross section for a few nuclei at about the same bombarding 
energy. As that determination was of a rather exploratory nature a remeasurement 
with improved accuracy was considered to be well motivated. 


2. Description of the experiment 
2.1. General considerations 


As in all attenuation experiments the equipment was designed to measure the 
number NV, of nearly monoenergetic protons incident on an attenuator (= target) 
and the number NV, of protons that emerge from the target without having caused 
any reaction. The total reaction cross section, o, can then be obtained from 


N,/Nq=exp (—n«) (1) 
N,—N 

‘tng for no<l, 
Ny 


or 
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The elastically scattered protons were distinguished from those inelastically scat- 
and from other reaction products by range separation. An absorber was inserted 
een the target and the detector for emerging particles. The energy resolution 
' hen necessarily limited by the combined effects of multiple Coulomb scattering 
and range straggling in this absorber. Target thicknesses and the geometrical proper- 
ies of the proton beam were chosen in such a way that their combined effect on 
nergy resolution should be small. Naturally some elastically scattered protons 
aused reactions in the absorber; i.e. the number actually counted was not N, but 
V.=%:N,, where 7 is the chance for an elastic proton to penetrate the absorber. 
‘The coefficient 1 was determined experimentally in a “‘target-out’’ measurement, 
oe the incident beam energy was reduced by an amount corresponding to the 
pereet thickness. This gave 
-where the prime denotes target-out measurements. Equation (1) then becomes 
e 2 No 
. Ny Nz 


= exp (—n0). (4) 


‘The target thicknesses were chosen so that the protons lost 9 MeV by ionization. This 
corresponded to no approximately equal to 0.02 and as 7 was about 0.79 a statistical 
accuracy of 1 % in the cross section would require a total number of incident protons 
of 1.4 x 10’. This could be reached in about 30 hours effective counting time, the 
limitation being the electronically allowable counting rate. 


2.2. The proton beam 


The experimental lay-out is shown in Fig. 1. As the requirements on the proton 
beam properties were very stringent the ordinary external beam from the cyclotron 
had to be modified. Four collimators were used at different positions along the proton 
path, marked A; B, C and D in Fig. 1. The collimator at B was a 10 mm thick brass 
piece with a circular aperture of 3 mm diameter. Its main purpose was to reduce the 
geometrical cross section of the beam. The thickness was not sufficient to stop protons 
from going through the wall, but on the other hand edge scattermg was kept small. 
At A the outer part of the proton beam was clipped, where the protons still had one 
turn left inside the cyclotron. The reason for the clipping was that it eliminated one 
of the two beams which would otherwise have appeared at C’. The collimators C and D 
were sufficiently large to let the main beam, i.e. protons coming through the aperture 
at B, pass without touching the edges. When the beam arrived at C it also contained 
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Fig. 1. The general lay-out of the experiment. 


a smooth background of low-energy particles. The small percentage of this tail that 
could pass through the aperture at C was made to miss the opening at D by the use 
of a bending magnet. Photographic methods were used to determine collimator 
sizes, and, in particular, heavy overexposures of photographic paper at D and at the 
exit of the last focusing magnet indicated that the beam was free from low-energy 
ionizing particles. Neutrons which were produced at C and D by the low-energy tail 
were negligible due to the very small solid angle of the counter when seen from D. 
Also the neutrons produced at A and B and in the magnetic channels inside the 
cyclotron were negligible due to the large distance from the detecting apparatus and 
the shields used. 

The beam was focused on the target by two pairs of quadrupole magnets, both of 
which were necessary to handle the large astigmatism of the beam emerging from the 
cyclotron. At optimum focusing the beam cross section at the entrance of the appara- 
tus was approximately circular with a diameter of 0.7 mm. The maximum angular 
spread obtained from the picture taken at the exit of the last focusing magnet, 
where the beam was also almost circular and had a diameter of 6-7 mm, was + 0.09°. 
This was completely negligible in comparison with the rms angular spread due to 
multiple Coulomb scattering in the Mylar foil, in the pre-absorber and in the first 
two counters which added up to about 0.4°. The size of the beam at the target was, 
however, essentially unaffected by this spread because the counters were very close 
to the target. 

The beam was always carefully aligned to pass through the centre of the detecting 
apparatus. The photographic methods used allowed a determination of the beam 
position to at least +0.2 mm at the scatterer and of its angular deviation from the 
symmetry axis to + 0.04°. 

During the runs the beam intensity had to be reduced in order not to overload the 


electronic circuits. In fact, the ion source of the cyclotron was turned off com- 
pletely. 
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Fig. 2. The detector arrangement. 


2.3. Apparatus 


_ The detector arrangements are shown in Fig. 2. The letters 4 and C refer to absor- 
bers of various materials and to scintillators (counters) respectively. The beam left 
he vacuum tube through a thin Mylar foil as close as possible to the first absorber. 

1 and C2 in coincidence gave directly the number of incident protons N,. Their 
ensions were 10 x 10 x 3 mm3 and 7 x 7 x 3 mm} respectively. The 4 mm thick 
ounter O3 was circular with a diameter of 100 mm and its distance from the target 
a chosen that it covered scattering angles up to 25°. The number NV, was obtained 
as the number of coincidences between C2 and C3. The absorbers A2 and A3 were 
shaped to give the same separation level between elastic and inelastic scattering 
events for all scattering angles; i.e. at a separation level of AH MeV the absorber 
thickness, experienced by a particle which had been scattered at an arbitrary angle 


( < 25°) and lost AZ MeV to the target nucleus in excess of the recoil energy loss, was 


equal to its mean range. This condition could only be exactly satisfied for one point 
in the target which was, of course, taken to be the middle but because of the small 
dimension of the beam and the relatively thin targets it was approximately true for 
all possible places of scattering. 

The aluminium absorber 42 was 10 Me¥ thick in the forward direction and shaped 
to compensate for recoil energy loss whereas the main part of the energy reduction 
occurred in the copper absorber 43, which consequently was of the correct shape for 
an infinitely heavy target nucleus and a plane parallel A2. As the recoil losses are 
different for different nuclei, each target had to have a separate absorber A2; they 
were, in fact, mounted in pairs on a wheel with remote control for placing an arbitrary 
target in the beam. The absorber Al was variable in eight steps and made it possible 
to measure at various points on the range curve. It was made of polythene in order to 
minimize small-angle scattering from it. 
_ The light pulges from the scintillators (NE102) were guided onto RCA6810 A 
photomultipliers from which the pulses were fed into transistorized coincidence units 
and fast prescalers. The dead time of the circuits was sufficiently small (= 25 nsec) 
for two protons from consecutive r.f. bursts of the cyclotron to be resolved. The 
counting rate, however, had to be kept below about 200 protons/sec because of the 
finite chance of two protons to occur within the same r.f. burst, and consequently 
give a false ratio N,/N. The dependence of this ratio on counting rates of this order 
was tested but no effect was observed. The resolving times of the coincidence units 
were 9 and 15 nsec for Cl + C2 and C2 + 03 respectively, which was found quite 
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sufficient to eliminate all accidental coincidences. The constancy of the efficiency 
of C3 over the large scintillator area was checked by determining the ratio N. 2/ N 
when 03 was being displaced vertically and horizontally and with all three absorbe ‘ 
removed. The ratio was 0.999 + 0.001 at all positions of 03, all displacements tried 


being less than 35 mm. 


2.4. Experimental procedure 


A series of runs was always begun by adjusting the cyclotron magnetic field to 
such a value that the extracted protons were of the correct energy. As reference 
the ratio N,/N, was used for the carbon target at a point where the range curve 
was very steep. The mean energy of the protons could by this method quickly be 
read to about 0.01 MeV (or 1 part in 18000) in relative measure. After that the 
beam was centered on the target by adjustments.of the bending magnet field for side 
variations and the cyclotron magnetic median plane for height variations. The beam 
position was checked by exposing photographic paper. 

As described in section 2.1., the actual runs consisted of determinations of the ratio 
N,/N, for two situations, namely “target in” and “target out”. Measurements were 
made at eight different positions along the range curve, i.e. for eight different thick- 
nesses of Al. Three positions were chosen to be well above the knee; they constituted 
the actual measurements of the reaction cross section for different levels of inelastic 
separation. The next three positions were used to determine the slope of the range 
curve, which had to be very well known in order to correct for shape effects above 
the knee. The range curve levelled off before reaching zero as Al was increased. The 
last two positions were chosen to determine that level. In the “target-out’’ measure- 
ments Al was always increased by a 9 MeV thick polythene piece to compensate for 
the energy loss in the target. The coefficient 7 was thus measured at the proper 
energy of the protons penetrating A2 and A3. As the shape of A2 was slightly diffe- 
rent for different targets the determination was made for two target-out situations 
(the carbon and the gold cases) in order to see if the size of the proton beam affected 
the N2/No ratio. No effect was found within the statistical errors. 

Errors caused by drift in the apparatus were minimized by frequent changes of 
targets and Al thicknesses. As the mean proton energy had to be known to about 
0.01 MeV the time variation of the cyclotron magnetic field was not always negligible. 
The energy reference mentioned above was frequently used during the measurements 
and corrections were applied in the calculations. 


3. Evaluation of the cross sections 


3.1. The shape of the range curves 


The use of equation (4) requires that the target-in and target-out measurements 
are made at “‘corresponding points” on the two range curves, i.e. at such points that 
the two measurements would have been equal if no reactions had occurred in the 
target. The shapes of the range curves are expected to vary with the target used 
even in the ideal case of no nuclear reactions either in the target or in the absorbers. 
As the slopes at measuring points even well above the knee of the curve were not 
negligible due to the extreme accuracy needed (about 1 part in 104 in the ratio 
N,/N, for an uncertainty in the cross section of 1%), the different shapes had to be 
taken into account. 
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Experimentally, this was done by defining corresponding points as the points for 
ich the distances along the range axis to the 50 % level were inversely proportional 
the slopes at that level. A consistent use of corresponding points for comparison 
target-in and target-out measurements not only eliminated first order shape 
cts but also automatically corrected for small deviations of the target thicknesses 
om 9 MeV as well as slightly different proton energies in the various runs. 
__-A theoretical estimate of the errors that could be expected from such a procedure 
was made. The shape of a range curve near its end is mainly affected by spread in 
energy of the incident protons, range straggling, multiple small angle scattering and 
geometry. 
_ (1) The spread in energy of the incident protons was difficult to estimate, but should 
not be larger than 1 MeV, which is known to be the upper limit for the energy spread 
of the ordinary external beam, i.e. the beam without collimators at A, B and O in 
Fig. 1. The exact value was not very important as its contribution to the slope of the 
) pene curve was small. In the calculations a total width of 0.6 MeV at half height was 
used. 

- (2) The range straggling was assumed to follow a Gaussian distribution with a rms 
‘Spread in range AR (g-cm-?) given by the formula (see for instance [10]). 


(A R)? =0.1571 “{(-3) “aT, (5) 


where 7 and 7’, denote the energies in MeV of the incident and the emerging protons 
respectively, (—d7'/dé) is the energy loss in MeV/g-cm~? and Z and A refer to the 
-charge and mass numbers of the material traversed, The values for (—dT'/dé) 
were taken from the tables of Aron [11]. The range straggling effect was dominant 
‘and corresponded to an energy spread of the incident beam of about 2.7 MeV full 
width at half height. The contributions to (AR)? from the different materials along 
the proton path are simply additive and the range straggling effect was also mainly 
responsible for the differences between the slopes at the 50% level when various 
targets were used. 

(3) A large number of small-angle deflections occur along the paths of the protons 
and make the actual path lengths greater than the thickness of the absorber. The 
theory of multiple Coulomb scattering by Moliére [12] was used to estimate this 
effect, which was almost completely determined by the thickness of the copper 
absorber 43, as small variations of the scattering angle in the target, did not appre- 
ciably alter the path length through the absorbers because of their curved surfaces. 
The path length distribution caused by multiple scattering alone was calculated by 
the method of Bichsel and Uehling [13]. 

(4) Several geometrical effects might contribute to the slope of the range curve. 
For instance, thé energy difference of protons leaving the target after being scattered 
elastically through 25° at the beginning or at the end of the target was 0.9 MeV. 
Beam width, angular spread, and misalignment should also have some effect. How- 
ever, all geometrical contributions were expected to be small and regarded as in- 
cluded in the energy spread of the incident protons. 

Assuming no inelastic events to take place but including the effects of range stragg- 
ling and multiple Coulomb scattering, the theoretical range curves were calculated 
for the C and Au targets by the method described by Bichsel and Uehling [13]. In 
the range straggling parameter was also included the effect of the energy spread of 
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because onl portant for the choice of correspondin 
Both theoretically and experimentally the slope was 2% smaller in the 
case than with the C target. } LouTeiRa 

To test the correctness of the method to choose corresponding points, the: 
axis of the theoretical Au-curve was changed in such a way that the slope at 
height coincided with that of the C-curve. Some small but significant deviations wi 
found between the C-curve and the transformed Au-curve. These deviations were 
used to find the ‘corrections due to shape variations” listed in Table 2. 


©, 
3.2. False coincidences in N, 


It was found that the N,/N, ratio levelled off at about 5 x 10-% for increasing Al 
thicknesses instead of going down to zero. One explanation for this fact would be 
that some protons after passing through C1 and C2 undergo (p,n)-reactions, where 
the neutrons proceed in the forward direction and have an appreciable chance of 
being detected by C3 either by ionization in the scintillator itself or by producing 
protons which can penetrate to C3. A rough theoretical estimate indicated that this 
effect might be of the order observed. Possibly also some decay gammas from nuclei 
excited by the protons might be detected by C3 and thus contribute to N,. The effect 
was not due to accidental coincidences or beam admixtures of either neutrons or 
protons of exceptionally high energy. | 

Because of the flatness of the low level it was considered reasonable to assume 
that the effect produced a constant background in N, of coincidences not resulting 
from direct protons. The size of the effect was determined for each target separately 
and for the two target-out cases. The target-out background was slightly smaller than 
that with a target in position but no other significant differences were found. The 
number N, to be used in equation (4) was in all cases subtracted by the proper back- 
ground determined in this manner. 


3.3. Scattering corrections 


Counter C3 covered scattering angles up to 25° and consequently the results had to 
be corrected for elastic scattering at larger angles and for inelastic scattering within 
that angle with so small an energy loss that the range separation used provided 
insufficient energy resolution. The corrections for elastic scattering were easily 
obtained by integrating the differential cross sections published in the following 
paper [1]. The uncertainties are mainly due to the poor knowledge of the limiting 
angle. 

It was considerably more difficult to estimate how many inelastically scattered 
protons could reach C3. The corrections applied were based on published inelastic 
differential cross sections [14] for those elements for which data are available. In the 
other cases corrections were found by studying the spectra obtained in an experi- 
ment designed to measure the elastic cross sections [1]. 
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Fig. 3. The total reaction cross section for protons at bombarding energies between 30 and 900 
; : MeV. References are given in Table 1. 
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Table 2. 


All values given are in 


Correction due to | 
. 
Level of fi : 
inelastic Uncorrecte % : ‘ energy d 
paoncny separation | cross section inelastic a shape ence bety 
(MeV) neattening oe variations | correspon 
detected at O>2 point 
pee fe A ee ae NS a ee eae 
~2. = 0 
Li 4.87 141.5 3:64 1.8 1.82 0.5 0 
3.53 153.3 2.84 1.4 =1.82:0:5 0 og 
3.02 146.4 2:54 fs3 =1.820:5 0 0 
Be 4.50 183.4 4.04 2.0 —1.2+0.3 0 0 
3.19 182.3 3.2+ 1.6 age oii). 3 0 0 
2.67 184.3 2.5 id? = be 2eniO.s 0 0 
Cc 4.41 204.5 16x 0.4 —1440.4 0 0 
3.10 216.2 0.84 0.2 —142+0.4 0 0 
2.59 215.5 0.5 0.1 —1.420.4 0 0 
Al 4.48 371.1 4.64 2.3 = 2.60.7 0.44 0.2 0 
3.15 401.6 3.64 1.8 = 2.6 OF 0.82 0.4 0 
2.63 387.8. 2.64 1.3 = 2-626 O21 0.44 0.2 0 
Ca 4.47 511.3 3.2+ 0.8 — 2.30.6 14+ 0.7 0 
3.13 530.4 125053 —2.3+ 0.6 2.0725. 1.0 0 
2.61 522.9 0.9 0.2 = 2.9 2 O16 2.0 1.0 0 
Fe 4.29 600 Dione aee iG a ick heel 0 
2.96 686 8 t 4 = foes | Bo ote = 2 0 
2.44 660 GM=eT3 Spee Hy eit l 0 
In 4.12 1093 22 UE rises a | 10” E=5 4+2 
2.79 1143 14 4 7 a eel 12) = 6 2+] 
2.27 1198 LOM eb =o) eae D 10m y S 225 
Au 4.37 1604 ie ena | Si ete 24 +12 8t4 
3.04 1598 LORRESS Sth shal) 32 +16 5+3 
2.52 1668 8 ser4 a ease | 28 +14 442 


4. Results and discussion 


The final results for the proton total nuclear reaction cross sections are given in the 
last column of Table 2. The mean proton energy was 179.6 MeV as obtained from the 
range-energy tables of Aron [11]. 

Fig. 3 shows the measured values of the cross sections from different experiments 
as a function of the bombarding energy between 30 and 900 MeV (note the break 
in the energy axis). With the exception of the 34 MeV data the results indicate that 
the energy dependence is slight within this energy band. 

The reaction cross sections from this experiment are compared with Kirschbaum’s 
results at 185 MeV in Fig. 4. As the agreement is quite good it seems that the 20 MeV 
separation level for inelastic events set by Kirschbaum was compensated by his use 
of a third counter which was sensitive only up to 20° scattering angle. 

In Fig. 4 are also plotted the values for the reaction cross sections obtained in an 
optical model calculation. The parameters of this model were derived from the 
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hions and errors. 
‘specified otherwise. 
2 Errors (standard deviations) from 


Averaged cross 


certainty of uncertainty in yeeaeereon section and total 
low-level position of corre- | counting statistics rms errors 
ubtraction sponding points 
0.8 0.2 4.6 143.3 
0.8 0.9 3.7 154.3 149.0+ 3.0 
0.8 1.5 4.0 147.1 
1.3 0.3 5.6 186.2 
: 1.3 1.3 7.2 184.3 185.5+ 4.5 
1.3 2.3 7.4 185.6 
AB 0.6 8.5 204.7 
1.5 2.7 6.3 215.6 212.3+ 4.8 
1.5 4.8 6.6 214.7 
3.0 1.2 17.6 373.5 
3.0 5.1 13.9 403.4 390 +10 
3.0 8.9 14.2 388.2 
4.4 1.8 25.2 513.6 
4.4 8.0 IOs) 531.3 524 +14 
4.4 13.8 21.7 523.5 
oe 5. 3 31 612 
mes 5 11 24 695 662 +19 
5 21 27 667 
9 5 52 1126 
9 28 42 1168 1165 434 
9 43 45 1217 
14 7 81 1646 
14 28 65 1641 1660 +50 
14 52 68 ~ 1704 


Gammel-Thaler potential for nucleon-nucleon interactions along the lines described 
by Kerman, McManus and Thaler [15]. The nuclear form factors were obtained from 
high energy electron scattering data [16]. The derivation is considered in more detail 
in the following paper [1]. It is evident that this optical model predicts too large 
cross sections except possibly for the heaviest elements studied. An approximate 
calculation showed that a decreased absolute value of the imaginary part of the 
potential would’ considerably reduce the predicted reaction cross sections for the 
lightest elements whereas the heaviest elements would be almost unaffected as they 
are already nearly opaque. 

The effect of nuclear correlations on the imaginary part of the potential is expected 
to be proportional to (V22 — Var) (see for instance [15)). Vor and Vc; refer to the 
real and imaginary central parts of the optical potential. As Vcg and Vc; are almost 
equal at 180 MeV it is unlikely that the discrepancy between experimental and theo- 
retical reaction cross sections can be explained by the neglect of nuclear correlations 
in the derivation of the potential parameters. 
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Fig. 4. The total reaction cross section for protons at 180 MeV (this experiment) and 185 MeV 
(Kirschbaum [6]) compared with the results of an optical model calculation. A is the mass number 
of the bombarded nucleus. ; 


On the other hand it is well-known that the Gammel-Thaler scattering matrix 
does not fit the proton-proton scattering data extremely well at energies around 
180 MeV and that the experimental information on the proton-neutron interaction 
is meagre. By using a more exact nucleon-nucleon scattering matrix one might 
therefore hope to be able to remove the discrepancy found or, conversely, the values 


obtained for the reaction cross sections might be useful for the study of the nucleon- 
nucleon interaction. 
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